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The optimal path crack model on uncorrelated surfaces, recently introduced by Andrade et al. 
{Phys. Rev. Lett. 103, 225503, 2009), is studied in detail and its main percolation exponents 
computed. In addition to 13 /u = 0.46 it 0.03 we report, for the first time, j/u = 1.3 ± 0.2 and 
T = 2.3 =b 0.2. The analysis is extended to surfaces with spatial long-range power-law correlations, 
where non-universal fractal dimensions are obtained when the degree of correlation is varied. The 
model is also considered on a three-dimensional lattice, where the main crack is found to be a surface 
with a fractal dimension of 2.46 =b 0.05. 
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I. INTRODUCTION 

Finding the optimal path between two points in a dis- 
ordered system is a relevant challenge for science and 
technology pTFTT]. This optimization problem is present 
in our daily lives when, for example, we make use of the 
Global Positioning System (GPS) to trace the best route 
to arrive to our destination. This problem, however, is 
not only relevant for human transportation. In materials 
science, the characterization of the optimal path is of ex- 
treme importance to study fractures, polymers in random 
environments, and transport in porous media [12 . If in- 
tensively used, this path is prone to fail, and a new path 
needs to be found. Studying how the successive paths 
evolve until the final configuration - where connectivity 
is no longer possible - is a challenge in itself which we 
address in this work. 

Recently, Andrade et al. [13 introduced a new model, 
named optimal path crack (OPC), to study the evolu- 
tion of successive optimal paths under constant failure. 
They have shown that, if a disordered energy landscape 
is considered, and each optimal path fails at its maximum 
energy site, the cracking process leads to a configuration 
where no more paths can be found. For uncorrelated 
systems, regardless the degree of disorder, the shortest 
path in the main crack - the minimal one to interrupt 
the connectivity - is always a self-similar object with a 
fractal dimension of 1.22 ± 0.02. This fractal dimension 
has been reported in several different systems like, e.g., 
the watershed line [14] and the perimeter of the percola- 
tive cluster at a discontinuous transition [15]. In this 
work, we first study in detail the properties of the OPC 
in uncorrelated lattices. The fractal dimensions are ac- 
curately obtained and, being a percolation-like process. 
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the main critical exponents of percolation are computed 
for the crack at the final adsorbing state [I6l [17] . Be- 
sides the previously reported value of in this work 
we show that j/u = 1.3 ± 0.2 and the Fisher exponent 
isr = 2.3±0.2. We then calculate and discuss the main 
properties of the model in three dimensions. 

The interplay between correlation and randomness lies 
at the very core of emergent phenomena. Here we also 
extend the study of the OPC to correlated lattices. As in 
many previous studies [T8]-[27], spatial long-range corre- 
lated energy distributions on the lattice are accessed by 
fractional Brownian motion (FBM) [28] - a generalized 
version of the classical Brownian motion, introduced by 
Mandelbrot and Ness [29] - where the degree of correla- 
tion between the successive steps can be tuned. 

The manuscript is organized as follows. In Section [Tlj 
the optimal path crack model is revised and its extension 
to correlated lattices is introduced. The results for uncor- 
related lattices, in the weak and strong disorder regime. 



are presented and discussed in Sec. |III| together with the 
main results for the three-dimensional system. The be- 
havior of the optimal path crack in correlated lattices is 
analyzed in Sec. IV We leave the final remarks for Sec.[V| 



II. MODEL AND DEFINITIONS 
A. The Optimal Path 

The optimal path crack (OPC) is obtained by system- 
atically interrupting the optimal path (OP) on a random 
landscape. 

For simplicity, let us consider a simple square lattice 
of linear size L where an energy > is assigned to 
every site. The energy of a path in a system is the sum 
over the energy of all sites in the path. The optimal path 
is the one, among all paths connecting two sites - e.g., 
the bottom and the top of the lattice - with the lowest 
energy. In the limit where all sites have the same energy, 
the OP is a straight line with mass (number of sites) 
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FIG. 1. (Color online) Double logarithmic plot of the mass 
of the optimal path, Mop, in 2D as a function of the linear 
system size, L, for two different types of energy distribution. 
The lower points (•) are measured for uniformly random site 
energies, being equivalent to weak disorder. The correspond- 
ing exponent is dop = 1.01 ±0.01 and the optimal path is a 
self-affine object [30]. The upper line (■) shows results ob- 
tained for site energies distributed according to a power-law, 
Eq. ([T]), with /3d = 1600, such that the systems are in the 
strong disorder limit (for the considered system sizes). In the 
strong disorder limit the optimal path is a self-similar object 
with a fractal dimension dop = 1.215 ± 0.005. The optimal 
path search was performed using the Dijkstra algorithm [31]. 
Results have been averaged over 10^° samples for the smallest 
system size and 10^ for the largest one. The error bars of the 
individual points are smaller than the symbols. 

When the energies are randomly distributed, the OP 
is, in general, not a straight line - its mass is larger or 
equal to L - and its properties depend on the energy 
distribution [30, 32-34 . In this work we take two energy 
distributions: uniform and infinite disorder (power law 
with exponent —1). In the former case, the mass of the 
OP, Mop, scales linearly with the system size as seen in 
Fig.[l] while this is not the case in the infinite disorder 
limit. 

When the special case of a power-law distribution with 
exponent —1 is considered, 

p{ei) (X — , (1) 

truncated between £max = 1 and Smin = exp(— /^i:)), the 
disorder in the system is controlled by the disorder pa- 
rameter Pjj > 0. In the limit Pd 0, Smin ^max, 

and the OP approaches a straight line. For large values 
of /3d, the truncated distribution becomes very broad 
and the disorder strong (ultrametric limit). The total 
energy of a given path is then mainly driven by the site 
with the highest energy in the path (see |32j and ref- 
erences therein). For a certain /3d ^ the system can be 
either in a strong or weak disorder regime, depending 
on its size. The crossover from weak to strong disorder 



has been investigated in detail by Porto et al [30]. For 
strong disorder, the OP is a self-similar object with a 
fractal dimension dop = 1.215 ± 0.005 (see Fig.[T|. 

The optimal path is typically the one heavily used in 
a system being liable to fail. This failure is most likely 
to occur in the highest energy site. Under failure, this 
site is destroyed and can no longer be used. To optimize 
the transportation "cost" on the modified substrate, a 
new OP has to be found, which is again disrupted at its 
highest energy site. During this process, transport in the 
system remains possible, with increasing cost, until the 
formation of a path of destroyed sites which disconnects 
the system into two parts. This process, denoted as op- 
timal path cracking^ was introduced recently by Andrade 
et al [13 , who discovered that the shortest path of de- 
stroyed sites necessary to disconnect the system has a 
fractal dimension of 1.22 ± 0.02, an interesting exponent 
also found in several other systems [14| \TE\ \30\ \32 \ \33\ l35] . 



B. Optimal Path Cracking 

The starting configuration for each realization of OPC 
is a regular square or cubic lattice of linear size L where 
all sites and bonds are occupied (open). Notwithstand- 
ing the site definition of the problem, the model can be 
straightforwardly extended to consider bonds instead of 
sites. Periodic boundary conditions are applied in all 
directions, except for one, where fixed boundary condi- 
tions are taken. For definiteness, this direction is called 
vertical, such that there are no bonds directly connect- 
ing the lowest row (layer) and the top row (layer) in two 
(three) dimensions. For each realization, a non-negative 
energy value e^, i = 1, . . . , = (where d is the di- 
mension of the system), is assigned to each site, either 
distributed uniformly (e.g. in [0,1)) or according to a 
truncated power law (Eq. ([T])). In the latter case, apply- 
ing the transformation method for distributions yields 

ei=exp[/3D{xi-l)] , (2) 

where X'l IS a random number uniformly distributed in 
[0,1). 

Starting from this setup, the OPC is generated by the 
following procedure: 

1. Find the OP through the lattice, connecting the 
lowest row and the highest one (the lowest layer 
and the top layer in 3D). Only occupied sites can 
be part of a path and only nearest neighbors are 
considered to be connected. 

2. Identify and remove the most vulnerable site in the 
optimal path, i. e., the one with the highest energy. 

3. Repeat steps 1 and 2 until there is no path connect- 
ing the bottom row (layer) to the top row (layer). 
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4. The configuration of unoccupied sites obtained in 
this way is the OPC. 

It is convenient to use Hoshen-Kopelmann-like label- 
ing [36H38 to keep track of the properties of the clusters 
of removed sites, as the OPC emerges. The properties 
of the cracks are analyzed by monitoring three different 



• The number of all removed sites in the system, Mtot 
- the density p is then defined as p = Mtot/N; 

• The number of blocked sites forming the largest 
cluster disconnecting the system, which is called 
Mic; 

• The shortest path in the largest cluster, sufficient 
to disconnect the system, with a mass Mgp. 

In Fig. [2] we show snapshots of the OPC for uncorre- 
lated lattices with different degree of disorder. The same 
sequence of uniformly distributed numbers was used for 
all configurations. Although the distribution of cracks 
and the number of dead ends are significantly affected 
by the degree of disorder, the shortest path in the main 
crack does not change. For uncorrelated lattices, the 
OPC shortest path is always a self-similar object with 
the same fractal dimension regardless the degree of dis- 
order. 

To find the mass of the shortest path in the largest 
cluster, the burning method introduced by Herrmann et 
al. [39] is used. Since for the optimal path only the four 
nearest neighbors are considered to be connected and an 
OP consists of connected occupied sites, it is sufficient 
for the OPC to be connected by next-nearest neighbors. 
Therefore, for the investigation of the OPC properties, 
each site is considered to be directly connected to its 
four nearest and four next-nearest neighbors (in 3D six 
nearest and eight next-nearest neighbors). 

The three masses introduced above are expected to 
asymptotically scale with the linear system size, L, 

Mtot^L^''^' , M,,^L^'^ , Msp-L^- , (3) 
where 



d\c = d- 



(4) 



is the fractal dimension of the largest cluster, {3 is the 
critical exponent related to the order parameter, and v 
is the exponent related to the correlation length. In addi- 
tion, we measure the second moment of the OPC cluster 
size distribution n. 



M2 = Y1 



(5) 



where Ug is the number of clusters of size s per lattice 
site, and Sk is the size of cluster k. Percolation theory 
predicts that the second moment, excluding the contri- 
bution from the largest cluster, scales asymptotically as 
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where Smax is the size of the largest OPC cluster and 7 
is a critical exponent. For d not larger than the upper 
critical dimension, the hyperscaling relation. 



d = 



(7) 



holds. The scaling behavior of the OPC cluster size dis- 
tribution gives access to the Fisher exponent r. 



The scaling relation between /3, 7, and r reads [16] 



7 



(8) 



(9) 



Andrade et al {IT showed that the OPC depends on 
the value of the disorder parameter Pjj. Small values of 
e.g., (Sd < 1, lead to narrow energy distributions 
and thus weak disorder. 



III. UNCORRELATED LATTICES 

For uncorrelated OPC with uniform energy distribu- 
tion, the total mass of the OPC, Mtot, the mass of the 
largest cluster of blocked sites, Mic, and the mass of the 
shortest path, Mgp, as a function of the linear system 
size are shown in Fig.jsJ As proposed in Eq. ([3|, the 
functional relation of the masses with the system size is, 
within the error bars, given by power laws. The obtained 
exponents are cZtot = 2.00 ± 0.01, d\c = 1.54 ± 0.03, and 
(igp = 1.21 ± 0.02, consistent to the ones for = 0.002 
^13j . Thus, in the weak disorder regime, the three masses 
scale with different exponents. Furthermore, since the 
mass of all blocked sites scales linearly with the number 
of sites in the system, the density of cracks can be com- 
puted giving p = 0.10 ± 0.03. When, instead of the site 
with highest energy, the one with lowest energy is re- 
moved, the same results as for OPC in the weak disorder 
regime are obtained, independently of the value of fSn- 

The variation of the crossover with system size from 
weak to strong disorder can be observed in Fig.|4j through 
the behavior of the three masses as a function of the 
linear system size. Results are for the intermediate value 
of = 6. One can observe that for small system sizes 
the system is in the strong disorder regime and the three 
masses coincide. The crack is then localized in a single 
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FIG. 2. (Color online) Snapshots of optimal path crack clusters on uncorrelated lattices, generated with the same seed for 
the uniformly distributed pseudo-random numbers. Figures (a) and (b) show equivalent realizations of the weak disorder limit: 
in (a) the site energies are distributed according to a power-law, Eq. ([T]), with 13d — 0.002, while in (b) they are uniformly 
distributed. Figures (c) and (d) show snapshots for f^D — ^ and — 100, respectively, corresponding to the intermediate and 
the strong disorder regime. The shortest path in the largest optimal path crack cluster is shown in red (medium gray), the 
sites belonging to the largest cluster but not the shortest path are shown in blue (dark gray), and the sites belonging to other 
clusters than the largest one are shown in green (light gray) . Notice that the shortest path in the largest cluster is independent 
on the distribution and on the disorder strength. Furthermore, since for the optimal path search only nearest neighbors are 
considered to be connected, the optimal path crack clusters are connected via nearest and/or next-nearest neighbors. The 
optimal path crack clusters are trees by construction. For these pictures, the linear system size is 256. 
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FIG. 3. (Color online) Double logarithmic plot of the to- 
tal mass, Mtot (■), mass of the largest cluster, M\c (•), and 
mass of the shortest path, Msp (a), as a function of the lin- 
ear system size, L, for OPC with uniformly distributed site 
energies. This energy distribution is equivalent to the weak 
disorder limit. The three masses scale with the system size 
according to three different exponents: dtot = 2.0 ± 0.01, 



die 



1.54 ± 0.03, and 4i 



1.21 zb 0.02. For weak disor- 



der, the OPC consists of a largest cluster with dangling ends 
and many small isolated clusters (see Fig. 2(a) and (b)). Re- 
sults have been averaged over 10^ samples for the smallest 
system size and 10^ for the largest one. The error bars of the 
individual points are smaller than the symbols. 



line. For increasing L the three curves separate and tend 
towards their weak disorder behavior. Andrade et al. [13] 
showed that the system size where the crossover occurs, 
I/x, scales with the disorder parameter: ^ Pd^ ^ with 
a ~ 0.59. In contrast to the behavior of Mtot and Mic, 
the mass of the shortest path in the largest cluster is 
independent of the system disorder and scales with the 
same exponent of dgp = 1.21 ± 0.02 as for jS^ = 0.002 
and for uniform energy distribution. 

Three main features characterize the qualitative be- 
havior of OPC clusters in the weak disorder limit. Firstly, 
when searching the OP, only nearest neighbors are con- 
sidered to be connected and the OPC growth ends when 
there is no path connecting opposite sites of the system. 
Thus the OPC clusters are connected via nearest and 
next-nearest neighbors. Secondly, by construction, ev- 
ery site in an OPC cluster was, at some point, part of 
an OP, therefore there can be no extended loops which 
would "trap" sites not belonging to the OPC. Compact 
loops, i.e., four OPC sites at lattice positions (xc^o), 
(xo, + 1), (xo + 1, ^o), and (xq + 1, + 1), can, how- 
ever, arise in OPC clusters. An example for the tree- 

Fi- 



like structure of OPC clusters is shown in Fig.|2(a) 



nally, branching in OPC clusters is less frequent than, 
e.g., in loop- less standard or invasion percolation clus- 
ters [33l [34l [40l [41] . These properties cause the shortest 
path, in the largest OPC cluster, to be identical to the 
ordinary backbone in the strong disorder limit, where the 



FIG. 4. (Color online) Double logarithmic plot of the to- 
tal mass, Mtot (■), mass of the largest cluster, Mic (•), and 
mass of the shortest path, Msp (a), as a function of the lin- 
ear system size, L, for site energies distributed according to a 
power-law with disorder parameter /^d = 6. The behavior of 
the total mass and of the mass of the largest cluster reveals 
the crossover from strong disorder to weak disorder. For small 
system sizes - equivalent to high disorder - all three masses 
coincide, while for increasing system size - decreasing disor- 
der - three different exponents emerge, as indicated by the 
asymptotic dashed lines which have slopes of 2.00 (for Mtot) 
and 1.54 (for Mic). The solid straight line is a guide to the 
eye with a slope of d^p = 1.21 zb 0.02. For this intermediate 
value of the disorder there are less isolated clusters than for 
weak disorder and the large st clu ster seems to be dominated 
by its shortest path (see Fig. 2(c) ). The inset shows the same 
masses as the main plot, measured for a higher value of the 
disorder parameter: = 30. For almost the entire range 
of system sizes the three masses are identical, indicating that 
these systems are in the strong disorder regime (see Fig. 2(d) ). 
The slope of the line gives dsp = 1.215 ± 0.005. Results have 
been averaged over 10^ samples for the smallest system size 
and 10^ for the largest one. The error bars of the individual 
points are smaller than the symbols. 



backbone becomes localized (see Fig. |2(d)| ). 

Next we discuss further details for the weak and strong 
disorder limits. 



A. Weak Disorder 

Let us analyze the properties of uncorrelated OPC 
clusters in the weak disorder limit. The critical expo- 
nents are estimated using the data for uniform site energy 
distribution. 

In the previous section we obtained, using Eq. 
13 /v = 0.46 ± 0.03. According to Eq. (|6|, the ratio 
can be determined from the scaling behavior of the sec- 
ond moment of the cluster size distribution, excluding 
the contribution of the largest cluster; the corresponding 
plot is shown in the inset of Fig.[5j Asymptotically, the 
points follow a power law with exponent j/i^ = 1.3 dz 0.2. 
Within the error bars, the obtained results fulfill the hy- 
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FIG. 5. (Color online) Double logarithmic plot of the cluster 
size distribution and its second moment for energies uniformly 
distributed. The main figure shows Ns = X^^ll^ rir (s = 
1, 2, 4, . . . ) as a function of the cluster size, s, for linear system 
sizes L = 512 (■) and L = 1024 (•). The slope of the line 
gives 1— T = — 1.3di0.2. In the inset we see the second 
moment, M2 (a), of the cluster size distribution as a function 
of the system size, L. The slope of the line gives j/u = 
1.3 ± 0.2. Results have been averaged over 10^ samples for 
the smallest system size and 10^ for the largest one. The error 
bars of the individual points are smaller than the symbols. 

perscaling relation given by Eq. 

From the behavior of the cluster size distribution for 
large cluster sizes we can extract the Fisher exponent r, 
as shown in Fig.jSj The points can be fitted by a line 
with slope — 1.3 ± 0.2, such that r = 2.3 ± 0.2 (details in 
the caption). This is consistent with the obtained values 
for P/u and j/u and the scaling relation Eq. ([9|. 



B. Strong Disorder 

For strong disorder, in contrast to weak disorder, the 
OPC cluster consists, mainly, in a localized crack, such 
that Mtot ^ Mic ^ Msp (see Fig. |2(d)[ ). As an example, 
the mass scaling in OPC, with site energies distributed 
according to Eq. Q and /^i:) = 30, is shown in the inset of 
Fig.|4] For the smaller system sizes the three masses co- 
incide, while for increasing system sizes one can notice a 
crossover of the total mass, Mtot- Since Mtot grows more 
slowly for higher disorder (with exponent ^1.2 instead 
of ~ 2.0), the amount of OPs that have to be found to 
produce the OPC is smaller than in the weak disorder 
limit. Since the multiple OP searches correspond to the 
most time consuming part of the algorithm, the required 
computation effort tends to decrease with increasing dis- 
order. This in turn allows to determine the fractal dimen- 
sion of the shortest path in the largest cluster with higher 
precision; we have found d^p = 1.215 ± 0.005, which is, 
within the error bars, identical to the OP fractal dimen- 
sion (ion = 1.215 ±0.005. 
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FIG. 6. (Color online) Double logarithmic plot of the mass 
of the largest cluster, Mic (•), for uncorrelated 3D optimal 
path cracks, as a function of the linear system size, L. The 
slope of the line gives die = 2.46 it 0.05. A snapshot of one 
realization is shown in Fig. [71 Results have been averaged over 
10^ samples for the smallest system size and 50 for the largest 
one. The error bars of the individual points are smaller than 
the symbols. 




FIG. 7. (Color online) Snapshot of a representative configu- 
ration of the optimal path crack (OPC) in three dimensions, 
for a system with 128*^ lattice sites. Energies have been as- 
signed to each site in the system according to a truncated 
power-law with /3d = 100. In black are the sites where the 
OPC intercepts the borders of the box. 



We also consider the uncorrelated OPC in three di- 
mensions, on a simple cubic lattice. Geometrically, it 
is clear that while a 2D system can be disconnected by 
a line, a surface is needed in 3D. Therefore, in analogy 
to the 2D case, one expects die to lie between two and 
three. For site energies distributed according to Eq. ([T]) 
with = 100, the mass of the largest OPC cluster as 
a function of the linear system size is shown in Fig. [6] 
The obtained exponent is die = 2.46 ± 0.05. A snapshot 
of an OPC, in three dimensions, in the strong disorder 
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regime, is shown in Fig.[7| Alike the two-dimensional 
case, for three dimensions also two different regimes are 
observed. While for the weak disorder regime several 
isolated clusters are found, in the strong disorder one 
the cracks become localized in a single surface. How- 
ever, computational limitations solely allow to accurately 
study the latter case. Simulations of systems in the weak 
disorder regime are actually limited to small system sizes 
and few samples. 



IV. CORRELATED LATTICES 

The distributions of energies in the system have been, 
so far, considered spatially uncorrelated. In this sec- 
tion, we introduce a generalization of the OPC study, 
described above, to correlated energy landscapes. These 
landscapes have been obtained from fractional Brownian 
motion (FBM) [ 19, i20» ^ 29 . We show that spatial 
long-range correlations lead to non-universal fractal di- 
mensions. 



A. Fractional Brownian Motion by Spectral 
Synthesis 

An ordinary Brownian motion [42 is a stochastic pro- 
cess B{t) with the properties 



and 



{B{t)-B{t'))=0 



{iBit)-Bit')f)cx\t-t'\ 



(10) 



(11) 



where t denotes time, and increments are mutually inde- 
pendent if, and only if, their time intervals do not overlap. 
Given a parameter called Hurst exponent, a fractional 
Brownian motion, is a moving average of dB(t)^ 

in which past increments of B{t) are weighted by the 
kernel (t — [29]. This leads to the properties 



{BH{t)-BH{t'))=0 



and 



{{Bnit) - BH{t')f) ^\t-t 



i\2H 



(12) 



(13) 



For H = 1/2 ordinary Brownian motion is recovered. 
In contrast to ordinary Brownian motion, in FBM, for 
H ^ 1/2, the increments are correlated and the range 
of this correlation is infinite. It can be shown that the 
correlation between two increments of a FBM, Bnit)^ is 
positive if 1/2 < < 1 and negative if < < 1/2. 

To use the properties of FBMs to obtain a correlated 
energy landscape, it is convenient to apply the Fourier 
filtering method (FFM), based on the spectral synthesis 



[T9[ [201 128] . Ill a nutshell, the idea is to generate random 
Fourier coefficients, distributed according to a given den- 
sity, and to subsequently apply an inverse Fourier trans- 
form to obtain the energy landscape in the spatial do- 
main. It is known that, in one dimension, 1// noise - a 
process with spectral density S{f) oc 1//^^ - is equiva- 
lent to FBM with H = {pc - l)/2. Furthermore, it can 
also be shown that the spectral density of an FBM in d 
dimensions can be written as 



S(!i,...,h) 



-{2H+d) 



(14) 



such that = 2H^d. Therefore, in 2D, for 2 < < 3 
the increments are anti-correlated, /3c = 3 corresponds 
to = 1/2 (ordinary Brownian motion), and 3 < Pc < 
4 leads to positively correlated increments. For white 
noise, the spectral density S is constant, i.e., /3c = or 
H = — (i/2, and both the increments and the obtained 
lattices are uncorrelated. 

In practical terms, we obtain correlated site energies 
in the following way [28 . For the spectral synthesis, 
Fourier coefficients corresponding to the spectral den- 



sity of Eq. (14) are needed. The given spectrum trans- 
lates into conditions for the expectation of the absolute 
values of the Fourier coefficients. We generate these 
Fourier coefficients in reciprocal space through one ran- 
dom phase (in [0,27r)) and one amplitude distributed 
according to a norma l distribution and with maximum 
amplitude {\/kl + /c^)"^^, where ki are the frequency 
indices of the discrete Fourier transform. After applying 
the inverse Fourier transform, we have to normalize in the 
range from zero to one the spatial domain distribution to 
represent the correlated topology. Separately, we gener- 
ate several samples and compute the average variance a. 
Afterwards, for each sample, we truncate its correlated 
distribution such that: if the distribution value is smaller 
(larger) than — 3<j (Sa), we assign this energy value equal 
to —3(7 {3a). This truncation affects less than 0.3% of 
the distribution. 



B. Correlated Optimal Path Cracking 

In this section, we study the properties of OPC on 
correlated lattices. Energies distributed equivalent ly to 
the spectral density in Eq. ( [l4| ) follow a FBM where Pc 
is the spectral exponent. Figure [8] shows the density of 
removed sites, p = Mt^t/N^ as a function of the correla- 
tion exponent Pc The number of sites that need to fail 
to break the global connectivity in the system decreases 
with the degree of correlation. The stronger the correla- 
tions, the lower the density of removed sites. While in 
the uncorrelated case {/3c = 0), the density of removed 
sites shows no significant finite-size effects, when spatial 
correlations are taken into account, the larger the system 
size, the lower the density of removed sites. 
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FIG. 8. (Color online) Density of all removed sites for cor- 
related OPC as a function of the correlation parameter, fie- 
For /3c = l, c?/2, the FBM corresponds to white noise, 

i.e., the uncorrelated case where the density is independent 
of the system size. Data for the following system sizes are 
shown: L = 32 (■), L = 64 (•), L = 128 (a), L = 256 (t), 
L = 512 (♦), and L = 1024 (•). Results have been averaged 
over 3.2 x 10"^ samples for the smallest system size and 10^ 
for the largest one. The error bars of the individual points 
are smaller than the symbols. 




FIG. 9. (Color online) The fractal dimensions dtot (B), die (•), 

and dsp (a) as a function of the correlation parameter, fie- 
The plot shows that the long-range spatial correlation changes 
the exponents of the system. Results have been averaged over 
3.2 X 10^ samples for the smallest system size and 10^ for the 
largest one. The error bars of the individual points are smaller 
than the symbols. 



The dependence of the exponents, dtot, dic^ and dgp, 
on is shown in Fig.[9| This plot shows that the crit- 
ical exponents of correlated OPC are non-universal and 
change with the correlation parameter Pc- In the ab- 
sence of correlations, the values of the exponents dis- 
cussed before are recovered. For /3c > 0, the exponents 
decrease monotonically with the correlation parameter. 
No significant differences are found between correlated 



and anti-correlated regimes. In the limit of strong long- 
range correlations, for sufficiently high values of f3c the 
exponents seem to converge towards unity. 

The density of all removed sites decreases monotoni- 
cally with Pc and the exponents for the different masses 
converge to unity, for large values of Pc Therefore, the 
OPC becomes localized in a single line. 



In Fig. [To] we see snapshots of the OPC shortest path 
and the first OP in four different energy landscapes, 
where the third dimension is the energy of the site. The 
snapshot (a) is for the uncorrelated case. Both the OP 
and the OPC have many upward and downward tilts. In- 
creasing the correlation, (a)-(d), the OP tends to cross 
the valleys and the OPC tends to go mainly through the 
mountains. Even for the anti-correlated case (b), though 
a rough landscape is also obtained, mountains and valleys 
can be observed. 



V. FINAL REMARKS 

In this work we have studied in detail the optimal path 
crack (OPC), recently introduced by Andrade et al [13], 
for the weak and the strong disorder limits. An accurate 
value for the fractal dimension of the shortest path in the 
crack was obtained. We also computed, for the first time, 
the set of critical exponents, j3/v^^/v^ and r, to get more 
insight on the crack percolation properties, revealing an 
interesting set of exponents. We note that, within the er- 
ror bars, the exponents P/v and ^/v are consistent with 
the ones for the parallel direction of directed percolation 
[43] . The Fisher exponent, r, is compatible with the one 
for the cluster size distribution of subcritical invasion per- 
colation clusters between two sites [44] and the exponent 
of the distribution of areas in perturbed watersheds [14J. 
The OPC has also been analyzed for a three-dimensional 
system where the connectivity between opposite borders 
of the system is broken by a surface with a fractal dimen- 
sion of 2.46 ± 0.05. 

We generalized the OPC on correlated energy land- 
scapes generated by fractional Brownian motion with dif- 
ferent values of the correlation parameter Pc For differ- 
ent correlations non-universal exponents are obtained for 
the fractal dimension of the total mass of the crack, the 
size of the main crack, and the length of the shortest 
path. Moreover, a monotonic decrease of the exponents 
with the correlation parameter is observed together with 
a strong tendency to localize the crack in its shortest 
path. An interesting extension of this work would be to 
consider, for example, in the three-dimensional case, the 
role of correlations and to compute ^/v and r. The in- 
terplay between degree of disorder and the presence of 
long-range spatial correlations is still an open question. 
Besides, since the obtained fractal dimension of the OP 
and of the shortest path in the largest OPC cluster was 
found in several different models, it would be interesting 
to understand what would be the meaning oi^jv and r 
for these models. 



9 





(c) 



(d) 



FIG. 10. (Color online) The shortest path in the optimal path crack for correlated and uncorrelated energy landscapes. At 
each lattice site, given by two coordinates {x,y), the third coordinate is the energy at this point. Red regions (mountains), 
have energy close to unity, while blue regions (valleys) have energy close to zero. The shortest path in the optimal path crack 
is shown in dark blue (dark gray, left to right) and the first optimal path in dark green (medium gray, front to back). Four 
different correlation parameters, /3c , have been considered: a) /3c = (H = —d/2) corresponding to the uncorrelated case - 
white noise; b) /3c — 2.2 (anti-correlated) ; c) /3c = 3 {H —1/2) recovering the ordinary Brownian motion where the increments 
(but not the lattice) are uncorrelated; d) /3c = 3.8, corresponding to strongly correlated lattices. While the crack passes mainly 
through the mountain tops, the first optimal path passes mainly through the valleys - hidden in the noisy configurations (a) 
and (b). For these pictures, the linear system size is 1024. 
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